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Abstract 

Minimum  norm  quadratic  unbiased  estimation,  or  MINQUE,  is  a  method  for  improving  the 
variance  estimates  of  noise  components  in  a  Gauss-Markov  least-squares  problem.  This  study 
treats  a  simple  special  case:  estimating  the  two  noise  levels  of  a  clock  whose  phase  noise  is  the 
sum  of  white  FM  and  random  walk  FM.  Given  prior  estimates  of  the  noise  levels,  perhaps  from 
an  Allan  deviation  plot,  MINQUE  calculates  new  estimates  and  their  uncertainties.  Although 
the  original  MINQUE  calculation  on  N  data  takes  0(N2)  space  and  0(N3)  time,  it  can  be  done 
sequentially  in  bounded  space  and  O(N)  time.  The  method  is  applied  to  data  from  a  simulation 
and  from  a  comparison  of  two  hydrogen  masers. 


INTRODUCTION 

Having  observed  time  deviations  (or  “phase”)  x(t )  of  a  clock  or  pair  of  clocks,  one  may  wish  to  fit  the 

noise  levels  ha  of  a  power- law  spectral  model,  Sv  (/)  =  hafa  ,  where  y  =  dx/dt .  This  is  often 

done  graphically  by  identifying  regions  of  constant  slope  on  log-log  plots  of  frequency  stability 
deviations,  including  Allan  deviation,  Hadamard  deviation,  and  their  modified  forms  that  use  interval 
averages  x  instead  of  point  values  x(t).  Objective  statistical  procedures  for  estimating  the  ha  and  their 

uncertainties  have  been  devised.  Tryon  and  Jones  [1]  estimated  the  noise  levels  of  an  ensemble  of  several 
cesium  clocks  by  maximizing  the  likelihood  function  of  the  residuals  of  a  Kalman  filter  operating  on  the 
clock  differences.  In  the  “multivariance”  methods,  the  data  consist  of  estimates  of  several  different 
frequency  stability  variances  at  a  set  of  averaging  times:  Vemotte  et  al.  [2]  estimated  the  noise  levels  by 
a  weighted  least-squares  technique;  Walter  [3]  estimated  the  noise  levels  and  mean  frequency  drift  rate  by 
minimizing  a  cost  function  composed  of  certain  chi-squared  statistics. 

The  present  paper  is  an  initial  effort  to  exploit  a  general  method  of  noise  level  estimation,  called 
minimum  norm  quadratic  unbiased  estimation,  or  MINQUE,  invented  by  the  statistician  C.  R.  Rao  [4]. 
This  is  a  method  for  improving  the  estimates  of  noise  levels  in  a  linear  model  z  =  Xb  +  £ ,  where  z  is  a 
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vector  of  observations,  X  is  a  known  matrix,  b  is  a  an  unknown  constant  vector,  and  the  noise  vector  £ 

»  ^  it  •  •  • 

has  the  form  2_,  <T,Z; ,  where  cr  are  unknown  nonnegative  coefficients  and  the  Z;  are  independent 
noise  vectors  with  zero  means  and  known  covariance  matrices.  MINQUE  accepts  prior  estimates  of  <7, 
and  produces  new  estimates  by  a  calculation  that  is  invariant  to  b . 

In  the  present  application  the  term  Xb  is  absent;  the  noise  is  all  there  is.  We  treat  only  a  simple  two- 
parameter  Gaussian  noise  model,  in  particular  the  phase  of  a  clock  with  unknown  levels  of  independent 
white  FM  (WFM)  and  random  walk  FM  (RWFM)  noises,  sampled  with  period  r0 .  The  mean  drift  rate  is 
assumed  to  be  zero.  The  data  for  MINQUE  consist  of  the  second  increments  of  the  phase  samples,  along 
with  prior  estimates  of  the  noise  levels.  The  first  step  of  MINQUE  is  to  normalize  the  unknown  noise 
levels  by  the  prior  levels.  The  second  step  is  to  prewhiten  the  data  according  to  the  prior  noise  model. 
The  last  step  is  to  apply  the  MINQUE  criterion  to  calculate  new  estimates  of  the  noise  levels  (Theorem  1) 
and  their  uncertainties  (Theorem  2).  We  show  examples  using  data  from  simulations  and  hydrogen  maser 
measurements,  and  we  see  the  effect  of  changing  the  prior  noise  levels.  We  also  try  feeding  the  estimates 
back  into  MINQUE  repeatedly  as  new  priors.  The  iterated  estimates  and  their  uncertainties  appear  to 
converge  to  reasonable  limiting  values  independent  of  the  original  priors.  The  Appendix  shows  results 
from  a  Monte  Carlo  trial  of  the  iterated  estimates. 

The  original  MINQUE  calculation  is  a  batch  operation  that  does  arithmetic  on  NxN  matrices,  where  N 
is  the  number  of  data.  After  it  was  invented,  statisticians  discovered  shortcuts  [5-7]  for  carrying  out  the 
calculation  in  much  less  time  and  space.  For  the  case  treated  here,  the  entire  calculation  can  be  carried 
out  sequentially  in  0(  N)  time  and  little  space.  A  short  summary  of  the  algorithm  is  given  below. 


SETUP  FOR  A  TWO-PARAMETER  MINQUE 

The  data  are  modeled  as  a  stationary  process 

z(n)  =  GlZl(n)  +  <J2Z2(n),  n  =  l,...,N.  (1) 

where  cr  >0.  The  base  processes  Zi ( n )  are  assumed  to  be  independent  Gaussian  stationary  moving- 
average  (MA)  processes  of  form 

Zi{n)  =  diQvi(n)  +  envi  (n-\),  *  =  1,2,  (2) 

where  the  are  known  coefficients  and  the  vt  ( n )  are  independent  Gaussian  random  variables  with 
mean  zero  and  variance  one. 

Special  Case:  White  FM  +  Random  Walk  FM 

Let  a  clock  phase  be  modeled  as  x(t)  =  xx(t)  +  x2(t) ,  the  sum  of  WFM  and  RWFM  noises.  The 
corresponding  frequency  spectrum  is  Sv(f)  =  h0  +  h_2f~2.  We  want  to  estimate  h0  and  h_2  from  N  +  2 
samples  x(nT0) .  Let  z(n)  be  their  second  increments, 

z(n)  =  x(«r0)-2x((«  +  l)r0 )  +  x((n  +  2)r0),  n  =  \,...,N,  (3) 
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which  form  a  stationary  process.  Unsealed  dimensionless  versions  of  the  two  components  of  z(n)  can  be 
represented  by  the  MA  processes 


Zl(n)  =  vl(n)-vl(n- 1)  (WFM) 
Z2(n)  =  v2(n)  +  j3v2(n-\ )  (RWFM) 


where  [5  =  2  -  V3  .  The  correlation  of  Z,  ( n )  and  Z2  (n  - 1)  is  1/4;  this  is  the  crucial  property  of  Z,  that 

makes  it  equivalent  to  the  second  increment  of  an  integrated  continuous-time  random  walk.  Now  (1) 
holds,  where 


=  hh. 

0  2’ 


07  =  h_2 


3(1  +  P2)' 


(4) 


Normalizing  by  Prior  Estimates 

Returning  to  the  general  case,  we  want  to  estimate  (7,  and  07  from  the  data  z  ( 1 ) , . . . ,  z  ( /V ) .  The 
M1NQUE  method  requires  prior  estimates  5]  >  0 ,  which  are  used  to  normalize  the  unknown  a,  .  Let 


Z=?s  z,(n)  =  d-,Z,(«),  an  =  afiir 

Instead  of  (1)  and  (2),  we  work  with  the  assumptions 

z(«)  =  TiA(”)  +  Z2z2(”)^ 
zi(n)  =  ai0vl(n)  +  anv1(n-l). 


(5) 


(6) 

(7) 


The  aim  is  to  calculate  estimates  yf  of  yf  from  the  data  z(/;j . 
negative.)  Then  erf  is  estimated  by 


(It  is  possible  for  ff  to  come  out 


(8) 


In  order  to  estimate  and  correct  the  overall  scaling  of  the  prior  values  o  ,  we  give  both  yi  an  unknown 
prior  value  £  that  will  be  estimated  from  the  data.  Expectation  and  variance  under  the  prior  assumptions 
7i~  Yi  =C  are  denoted  by  Ep  and  varp . 


Prewhitening  the  Data 

The  data  z(«)  will  be  transformed  into  new  data  y(n)  that  are  white  under  the  prior  assumptions.  Write 

z,.=[z,.(l)  ...  z,.(A)]T,  i  =  1,2,  (9) 

Z=7lZl+72Z2‘  (10) 


Define  the  Nx(N  +  \)  Toeplitz  matrix 
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L  = 


0-n  CXio 


(ID 


and  let  T  =  LXL]  +  L2Lt2  ,  which  is  assumed  to  be  positive  definite.  By  (7),  (10),  and  the  prior  assumptions, 

E  zizJ=LiUi,  EpzzT=C2T. 

By  Cholesky  factorization,  T  =  LL'  for  a  nonsingular  lower-triangular  L,  which  we  use  to  reduce  z  to  an 
approximation  of  white  noise.  Let  y  =  P'z  ,  v,  =  P'zi .  From  (10), 


y  =  YP\  +  YiYi » 


and  yx,y2  are  independent  mean-zero  Gaussian  vectors.  Let 


Then 

where 

The  true  covariance  matrix  of y  is 


Vi=EyiyJ. 

Vi  =  L-'L^P1  =MiMj  , 
M,  =  zTz., . 

EyyT  =Y\V\  +  YiV2  ■ 


But,  since  Vx  +  V2  =  L  lTL  1  =  / ,  y  is  white  under  the  prior  assumptions: 


>yy 


=  C2i- 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 


MINQUE  ESTIMATORS 

Matrix  Preliminaries 

1.  Let  A  and  B  be  real  matrices  of  size  m  x  n  .  Regarded  as  vectors  with  mn  components,  their  inner 
product  is 

{A,B)  =  YJajbij  =  trABT, 


where  tr  is  the  trace  operator.  The  Frobenius  norm  of  a  matrix  A  is 

\\A\\  =  yl(A,A)  =  =sltrAAT  . 
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2.  If  A  is  symmetric  and  y  is  a  column  vector  of  correlated  Gaussian  random  variables  with  E y  =  0  , 
E  yyT  =  V  =  [v,y] ,  then 


E yT Ay  =  (A, V)  =  tr  AV ,  (18) 

war  yT  Ay  =  2  (AV,VA)  =  2trAVAV.  (19) 

Equation  (19)  can  be  proved  using  Isserlis’s  formula  [8], 

cov(y[yJ,ykyt )  =  vikvjt  +  vuvjk . 


MINQUE  Criterion 

A  MINQUE  estimator  of  yf  (i  =  1,2)  from  the  prewhitened  data  vector  y  of  (12)  is  defined  as  a  random 
variable  yf  =yTAiy,  where  Ai  is  a  symmetric  matrix  chosen  such  that  yf  is  unbiased  for  yf  and  has 
minimum  variance  under  the  prior  assumptions  yl  =  y2  =  f . 

MINQUE  Solution 

We  shall  see  that  the  MINQUE  estimators  do  not  depend  on  f  .  One  needs  to  know  only  the  data  vector 
v  and  the  covariance  matrices  Vi  =EyiyJ  that  come  from  the  prewhitening  calculation. 


THEOREM  1.  Assume  that  Vl,V2  are  linearly  independent  as  N  x  /V  matrices.  Let 


~(VX,VX)  {vx,v2) 

-yTKi 

_{V2,VX)  {v2,v2)_ 

>  I  - 

yv2y_ 

The  unique  MINQUE  estimators  of  yf ,  y%  are  given  by 


f:= 


ff 

n 


=  S~'q. 


(20) 


(21) 


Proof.  Let  Ax  be  symmetric.  By  (18)  with  A  =  At  and  V given  by  (16), 

E yTA,y  =  tf  (A, ,  Vx )  +  f2  (4 ,  V2 ) . 


By  (19)  with  A  =  Ax  and  V  given  by  (17)  under  the  prior  assumptions, 

varp  yJ  Axy  =  2f4\\A1f  .  (22) 

For  yT  Axy  to  be  unbiased  for  yf  whatever  the  true  values  of  yx  and  y2 ,  it  is  necessary  and  sufficient  that 

(4,^  =  1,  (4,v2)= 0,  (23) 
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which  is  possible  because  of  the  linear  independence  assumption.  According  to  the  M1NQUE  criterion 
and  (22),  we  are  supposed  to  minimize  the  matrix  norm  of  Ax  while  keeping  its  matrix  inner  products 

with  Vx  and  V2  fixed  at  1  and  0.  The  minimizing  Ax  is  necessarily  a  matrix  (an  N2  -vector)  in  the  span  of 
V}  and  V2 ;  any  component  orthogonal  to  that  span  would  increase  the  norm.  Thus,  Ax  =cnVx  +cnV2  for 
some  coefficients  cn,cl2.  A  similar  argument  for  yTA2y  gives  A2  =  c2]Vt  +c22V2  for  some  c21,c22 .  Let 
C  =  \_cy ]  •  Altogether  we  have  (^Ai,Vj^  =  S:j  (Kronecker  delta),  and  substituting  for  the  Al  gives  CS  =  I , 
C  =  S~l .  The  M1NQUE  estimates  of  the  yf  are  then 

K  =yT  Ay =yT(cnvi+ci2v2)y  =  cu‘h  +cu(h 

y2  —  y  A2y  —  y  (c21Fj  +  c22V2 ) y  —  c2xqx  +  c22q2. 


The  original  unknowns  of  are  estimated  by  (8)  from  the  prior  estimates  df  •  As  mentioned  earlier,  ff 
can  be  negative. 

Estimated  MINQUE  Covariance 

THEOREM  2.  Under  the  prior  assumptions  y,  =  y,  =  T  ,  the  MINQUE  estimator  y1  has  covariance  matrix 

R  =  2  C4S~l. 

Proof.  Let  Ai ,  S,  C  be  as  in  the  proof  of  Theorem  1 .  For  an  arbitrary  2-vector  u  =  \ux  u1  ]T  , 

u 1  y2  =  iqy1  Axy  +  u2yTA2y  =  yT  (ulAl  +  u2A2)y 
By  (17)  and  (19),  under  the  prior  assumptions, 

var vuYy2  =2f4jjulAl  +«2J2|‘  =2  f4'^uiuj(^Ai,AJ^. 

hj 

But 

(4.4)=Zv^^)=Iwi(  =(csc)  =c, 

k,l  k,l 

since  C  =  S~l .  Thus, 

varp«Ty2  =u 1  (2^45_1)« 

for  any  u.  It  follows  that  2^45_1  is  the  prior  covariance  matrix  of  y2 .  ■ 

Under  the  prior  assumptions,  yTy/C2  is  a  /f  random  variable,  so  let  us  estimate  f2  and  by 

C2=yTy/N,  Rf=2pS~K  (24) 

According  to  (8),  we  can  estimate  the  standard  deviation  of  &2  from  prior  assumptions  by 
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std  (25) 

Neither  the  estimates  (8)  of  of  nor  their  estimated  standard  deviations  (25)  change  when  the  prior 
estimates  ox  and  o2  are  multiplied  by  the  same  constant.  But  the  examples  given  below  show  that  (25) 
is  sensitive  to  multiplying  them  by  different  constants. 


EXAMPLES 

Simulated  Noise:  WFM  +  RWFM 

A  simulation  of  A  =  1000  values  of  z(n)  as  given  by  (3)  was  generated  with  the  parameters  r0  =  1, 
h0=  1 ,  h_ 2  =  1.900  xlO"4 .  Results  in  the  form  of  Allan  deviation, 


h_22n2z 

3 


(26) 


[9],  are  presented  in  Figures  1  and  2  for  different  prior  values  h0  and  h_2  mapped  to  erf  and  of  by  (4). 
The  black  curve  shows  the  Allan  deviation  for  the  true  ha .  The  figures  also  show  the  overlapped  Allan 

deviation  estimate  for  z  =  2"  with  68.3%  confidence  intervals  calculated  by  the  Stable  software  [10].  No 
confidence  interval  is  calculated  for  the  last  point. 

For  Figure  1  we  set  prior  estimates  h0  =h0/2,  h_2  =  2 h_2 ,  that  is,  yf  =  2 ,  yf  =1/2 .  These  give  the  prior 
Allan  deviation  shown  by  the  red  curve.  The  blue  curve  is  the  Allan  deviation  (26)  from  the  estimates  ha 
as  mapped  from  the  M1NQUE  estimates  of  by  (4).  The  dashed  blue  curves  are  the  Allan  deviations  for 
ha  +stdp  ha  and  ha—  stdp  ha ,  where  the  estimated  standard  deviation  stdp  ha  is  calculated  from  (25)  and 
(4).  The  lower  curve  terminates  when  it  goes  negative,  because  stdp  h_2  >  h_2 . 
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Figure  1.  MINQUE  estimate  of  simulated  WFM  +  RWFM  noise.  The  prior  noise  level 
was  set  to  half  the  true  level  for  WFM,  twice  for  RWFM.  The  error  bars  show  the 
overlapped  Allan  deviation  estimate. 


For  Figure  2,  we  set  h0  =  2 h0 ,  h_2  =  h_2  /  2 .  Now  the  ±1  stdp  limits  for  the  RWFM  component  seem  to  be 
too  tight.  For  both  cases,  the  WFM  estimate  hQ  and  its  standard  deviation  are  insensitive  to  the  choice  of 
prior  estimates.  The  RWFM  estimate  h_2  is  somewhat  more  sensitive  to  the  priors,  but  is  still  a  good 
estimate  of  h_2 .  Its  estimated  standard  deviation  is  very  sensitive  to  the  priors,  however. 


Figure  2.  The  ratios  of  prior  levels  to  true  levels  are  reversed. 
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Figure  3.  After  5  feedback  iterations  from  either  set  of  original  priors. 

Iterated  Feedback  of  Estimates 

Flaving  calculated  the  M1NQUE  estimates,  one  can  repeat  the  process  with  the  same  z  data  by  using  the 
ha  in  the  role  of  the  prior  ha  and  calculating  new  M1NQUE  estimates  of  ha .  This  process  can  be 
repeated  indefinitely  as  long  as  both  ha  are  positive.  For  the  examples  that  the  author  has  tried,  the 
estimated  ha  and  standard  deviations  appear  to  converge  to  limits  that  are  independent  of  the  original 
priors.  The  estimated  normalizing  factor  £  from  (24)  tends  quickly  to  1.  For  the  simulated  data,  Figure 
3  shows  the  result  after  5  iterations.  By  luck,  the  limiting  RWFM  estimate  is  close  to  the  true  value.  The 
standard  deviation  estimate  looks  reasonable  in  view  of  the  Allan  deviation  confidence  intervals.  There 
are  two  questions:  1)  Is  the  limiting  ha  unbiased?  2)  How  well  does  the  limiting  standard- deviation 
estimate  represent  the  standard  deviation  of  the  limiting  hal  To  answer  these  questions  empirically,  one 
can  collect  statistics  on  the  limiting  estimates  over  many  simulated  realizations  of  the  clock  noises  (see 
the  Appendix). 

Real  Phase  Noise:  Two  Hydrogen  Masers 

Figure  4  shows  overlapped  Allan  deviation  with  confidence  intervals,  the  MINQUE  estimate,  and  its 
estimated  ±  1  std  limits  for  a  2-week  comparison  of  two  hydrogen  masers  in  JPL’s  Frequency  Standards 
Laboratory.  The  phase  data  x(t)  were  subsampled  to  15  minutes  from  an  original  sample  period  of  1 
second  to  get  beyond  the  region  of  white  and  flicker  phase.  Their  second  increments  (3)  constitute  the 
data  for  MINQUE.  The  overlapped  Allan  deviation  estimate  with  error  bars  is  also  shown.  The  red  curve 
is  Allan  deviation  from  the  prior  noise  parameters  h0  =5xl0“27  s  ,  h_2  =7.45xl0“36s"1 ,  which  were  fit  by 
eye  from  the  Allan  deviation  estimates  of  the  original  data  set.  The  prior  WFM  level  seems  too  high  for 
the  subsampled  data,  and  MINQUE  adjusts  it  downward  appropriately.  Although  RWFM  looks  like  a 
poor  fit  for  the  long-term  behavior,  MINQUE  chooses  a  reasonable  estimate  of  the  noise  level.  After  five 
cycles  of  iterated  feedback,  the  estimates  do  not  change  enough  to  make  another  plot  worth  showing. 
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Figure  4.  Applying  M1NQUE  to  the  phase  difference  of  two  H-masers,  assuming  a  noise 
model  of  WFM  +  RWFM. 

SEQUENTIAL  MINQUE  ALGORITHM 

If  no  attention  is  paid  to  the  special  structure  of  the  matrices,  the  prewhitening  calculation  for  N  data  takes 
0{  N3 )  time  and  (9(/V2)  space.  The  author  was  able  to  devise  an  iteration  that  calculates  the  S  and  q 

matrices  of  Theorem  1  from  the  z  data  and  coefficients  aij  in  O(N)  time  and  bounded  space.  For 
77=1,2,...,  the  algorithm  receives  z(n)  and  updates  a  small  structure  of  scalars,  2-vectors,  and  2x2 
matrices  including  y(n),  S,  and  q.  The  details  are  available  from  the  author.  For  small  values  of  n, 
equation  (21)  can  give  negative  values  of  ff ;  this  is  of  no  concern  and  the  iteration  can  proceed. 


CONCLUSIONS 

Several  techniques  were  introduced  while  implementing  this  case  of  Rao’s  MINQUE: 

•  Applying  the  MINQUE  criterion  after  prewhitening  the  data  according  to  the  prior  estimates  of 
noise  levels; 

•  Introducing  and  estimating  a  scaling  factor  £  as  the  prior  value  of  the  normalized  noise  levels; 

•  Devising  an  O(N)  algorithm  for  the  calculation  (as  also  claimed  in  [6]); 

•  Estimating  the  noise  level  uncertainties; 

•  Feeding  the  estimated  noise  levels  back  into  the  algorithm  as  new  prior  levels. 

The  MINQUE  method  seems  to  work  well  for  estimating  the  noise  levels  of  the  simple  clock  model  of 
white  FM  plus  random  walk  FM  with  no  drift,  giving  lower  RWFM  uncertainty  than  overlapped  Allan 
deviation  does.  The  method  requires  prior  estimates  of  the  noise  levels  as  inputs;  for  best  results,  it  seems 
advisable  to  feed  the  output  estimates  back  into  the  algorithm  as  new  priors  at  least  once  or  twice.  The 
feedback  experiments  show  empirically  that  the  iteration  converges  to  a  fixed  point  not  dependent  on  the 
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original  priors.  The  fixed  point  constitutes  another  estimate  of  the  noise  levels.  Properties  of  this 
estimate  and  its  associated  uncertainty  matrix  should  be  investigated  theoretically  and  empirically  (see  the 
Appendix). 

The  sequential  version  of  the  algorithm  might  be  simplified  by  using  a  time-invariant  filter  for  the 
prewhitening  operation.  Perhaps  the  idea  of  applying  the  M1NQUE  criterion  to  prewhitened  data  could 
be  extended  to  more  complex  situations  with  several  noise  components  or  several  clocks. 
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APPENDIX:  EXPERIMENT  ON  ITERATED  MINQUE 

To  investigate  the  properties  of  the  fixed  point  of  the  MINQUE  feedback  iteration,  1000  independent 
realizations  of  1 000  values  of  WFM  +  RWFM  noise  were  generated  with  the  same  noise  levels  ha  as 
before.  For  each  realization,  five  feedback  MINQUE  iterations  were  run  from  randomized  initial  prior 
noise  levels  to  get  close  to  the  fixed  point.  Statistics  were  collected  on  the  noise  level  estimates  ha  and 
their  estimated  standard  deviations.  Figure  5  shows  the  results. 
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Figure  5.  Experimental  statistics  of  a  case  of  iterated  MINQUE. 


The  left-hand  plots  show  the  histograms  of  h0  (top  plot)  and  its  estimated  standard  deviation  (bottom 
plot),  normalized  by  the  true  value  h0 ;  the  right-hand  plots  show  similar  histograms  for  h_2 .  For  both 
noise  levels,  the  sample  mean  of  ha  is  close  to  the  true  value,  and  the  sample  standard  deviation  of  ha 
(std  value  in  top  plots)  is  close  to  the  sample  mean  of  the  standard  deviation  estimate  (mean  value  in 
bottom  plots).  The  standard  deviation  estimate  for  h_2  has  a  considerable  spread,  however. 

Rao  mentioned  the  possibility  of  MINQUE  iteration  [11]: 

These  estimates  may  then  be  substituted  in  (3.15),  and  the  MINQUE  procedure  repeated.  The  whole 
process  may  be  repeated  several  times,  but  then  the  property  of  unbiasedness  is  usually  lost.  It  is 
possible  that  the  estimates  so  obtained  may  have  other  interesting  properties. 

For  the  case  investigated  here,  the  iterated  noise  level  estimates  and  their  associated  standard  deviation 
estimates  appear  to  be  unbiased  for  practical  purposes. 
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